Time-of-flight observables and the formation of 
Mott domains of fermions and bosons on optical lattices 
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We study, using quantum Monte Carlo simulations, the energetics of the formation of Mott 
domains of fermions and bosons trapped on one-dimensional lattices. We show that, in both cases, 
the sum of kinetic and interaction energies exhibits minima when Mott domains appear in the trap. 
In addition, we examine the derivatives of the kinetic and interaction energies, and of their sum, 
which display clear signatures of the Mott transition. We discuss the relevance of these findings 
to time-of-flight experiments that could allow the detection of the metal-Mott-insulator transition 
in confined fermions on optical lattices, and support established results on the superfluid-Mott- 
insulator transition in confined bosons on optical lattices. 
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Following the successful loading of ultracold Bose gases 
on optical latticesji*^ a similar scenario has been achieved 
in recent experiments with ultracold fermions^*^*^ In 
Refs.|3 and0 an optical lattice was superposed on highly 
anisotropic traps to study transport properties and lo- 
calization effects in ideal Fermi gases and Bose-Fermi 
mixtures. More recently, Hubbard- and Luttinger-type 
systems have been realized by loading fermionic atoms 
on three— and two-dimensional^ optical lattices, respec- 
tively. Now that such experimental setups are available, 
one of the most prominent goals is the achievement of the 
metal-Mott-insulator transition^ in confined fermions on 
optical lattices&S 

Ultracold fermions on optical lattices are an almost 
ideal realization of the Hubbard model, in which all pa- 
rameters can be controlled with very high precision. In 
this work, we are interested in the one-dimensional (ID) 
regimefSiS where the Hamiltonian can be written as 
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where , Cio- are the creation and annihilation operators 
of a fermion in a pseudospin state a at site i (and posi- 
tion a;,;), and rii = J^a^icr {riia = cl^c^^) is the particle 
number operator. The on-site interaction parameter is 
denoted hy U {U > 0), the hopping amplitude by t, and 
V2 is the curvature of the harmonic confining potential. 
In experiments, t and U can be modified by changing the 
intensity of the laser beams that produce the lattice. In 
addition, U can be also controlled separately by means 
of a Feshbach resonance which is the mechanism we 
consider below. Our quantum Monte Carlo (QMC) sim- 
ulations of the fermionic Hubbard model were performed 
using a projector algorithmiSiii along the lines described 
in Refs.]^ and|^ 



In the absence of a confining potential (V2 = 0), the 
phase diagram of this model, Eq. for U > 0, consists 
of a Mott-insulating (MI) phase at half-filling (n = 1), 
trivial band insulating phases for n = 0, 2, and a metallic 
phase for any other density. In the presence of a con- 
fining potential V2 ^ 0, these phases coexist in space 
separated regionsi^ Local measurements may then pro- 
vide an experimental proof of the existence of MI do- 
mains in the trap^sS A recent proposaU^ also includes 
collective oscillations as a possible way to detect their 
formation. However, the demonstration of the MI for 
fermions is more difficult than in the bosonic case. This 
is because the appearance of Mott regions for bosons can 
be detected experimentally by means of standard time- 
of-fiight (TOF) measurements^*^ while in the fermionic 
case the momentum distribution does not exhibit any 
distinguishing feature when MI plateaus develop in the 
trapi^ 

We propose here a novel way to detect the metal-MI 
transition in fermionic systems. We find that the sum 
of kinetic and interaction energies, which can be mea- 
sured in TOF experiments, and their numerical deriva- 
tives show distinctive signatures of the formation of Mott 
domains of fermions as well as bosons trapped on optical 
lattices. Hence, energy measurements provide a unified 
way to deal with the Mott transition in confined atoms 
(fermions or bosons) on optical lattices. 

Bosons on ID optical lattices can be described by the 
boson Hubbard modeli^ 
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where a|, Ui are the creation and annihilation operators 
of a bosonic atom at site i, and rii = a|ai is the particle 
number operator. The hopping parameter, on-site inter- 
action, and curvature of the confining potential are de- 
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noted as in the fermionic case. The QMC simulations of 
the boson Hubbard model were done using the stochastic 
series expansion (SSE) methodic in the grand canonical 
ensemble in which an additional chemical potential term 
in Eq. (0) is used to control the filling. 

In Fig. n^a), we show the evolution of the variance of 
the density (A^ = (n^) — (rii)^) in a ID system of fermions 
when the on-site repulsion is increased using a Feshbach 
resonance. At J7 = 0, n,; < 2 throughout the trap so 
that no local "band insulator" is present. This can be 
easily checked experimentally measuring the momentum 
distribution function, which should exhibit a Fermi mo- 
mentum kp, i-e, Ukykp = Oii^ (The band insulator al- 
ways produces a finite occupation of all momentum states 
within the Brillouin zone. For U > At, a. region of small 
Ai, indicating suppressed (but finite) quantum fluctua- 
tions of the density, appears in the center of the trap 
showing the establishment of a MI thereiS*^ Notice that 
in a trap, in contrast to periodic systems at half-filling, 
the MI appears only for a finite value of U/t. 



in 



in 



□ 



10 



(a) 



0.1 0.2 0.3 0.4 0.5 



(b) 



-24 



-12 





x/a 



12 



24 -24 



-12 





x/a 



12 



24 



FIG. 1: (color online). Intensity plot of the density fluctua- 
tions Ai = {nj) — {rii)'^ vs position and interaction strength, 
in fermionic (a) and bosonic (b) systems with 30 atoms on 
48 lattice sites and V^a^ — 0.015t (a is the lattice constant). 
The lattice size was chosen large enough so that the particle 
density is zero at the boundaries. For large U, the constant 
small value of Ai in the center of the trap signals the forma- 
tion of the MI core. In (b), for low values of U/t, we have 
truncated A^ at 0.5, i.e., for low on-site repulsive interactions 
the bosonic fluctuations of the density exceed 0.5. 



at half maximumjSiiSiiS and more recently the behavior of 
the visibilityi^SiSi They are all related to the momentum 
distribution function Uk of the trapped system, which 
is measured in TOF experiments. We have studied these 
observables in the fermionic case, and found no signature 
of the formation of MI domains. 

A further quantity related to Uk that can be also ob- 
tained experimentally is the kinetic energy. 



Ek = — 2t^^n/c cos(A:a), 



(3) 



where a is the lattice parameter, and — 2icos(fca) is the 
lattice dispersion relation. 




FIG. 2: (color online), (a), (d) Kinetic energy (Ek), (b), (e) 
interaction energy (Ej), and (c), (f) the sum Ek + Ej, for 
fermions (a-c) and bosons (d-f ) in the systems of Fig. Q (30 
atoms and ¥20^ = 0.015t). Energies are given in units of t. 
Continuous lines depict numerical derivatives of the energies. 
The arrows signal the emergence of a MI core. The QMC 
errors are much smaller than our symbol sizes throughout 
this work. 



The corresponding evolution of A^, for a bosonic sys- 
tem with the same trap parameters and number of par- 
ticles as the fermionic case, Fig. ^a), is depicted in 
Fig. n^b). Ai again signals the formation of a MI 
plateauii when U ^ 6. It, i.e., a value that is ~ 2t larger 
than in the fermionic case and also larger than the crit- 
ical value U/t ~ 3.5 in the periodic casepi& Although 
for small values of U , quantum fluctuations of the den- 
sity for bosons are much larger than for fermions (in the 
figure we have truncated them at A = 0.5), when U is 
increased both systems behave analogously, revealing the 
similarities between bosons and fermions in ID. 

Several approaches based on TOF measurements have 
been followed in experiments to detect the formation of 
Mott domains of bosons. They include the disappearance 
of the interference pattern^ the increase of the full width 



The evolution of the kinetic energy with U, for the 
fermionic system in Fig.^ is shown in Fig.|2Ia). Ek does 
not exhibit any special feature at the point where the MI 
appears in the system. In contrast, the interaction energy 
(Ej), depicted in Fig. El^b), displays a sharp drop just 
before the Mott plateau sets in. After that, the decrease 
of Ej occurs slowly, even more slowly than the reduction 
of the modulus of the kinetic energy. This suggests that 
if one turns the trap off and lets the particles expand on 
the lattice, so that the interaction energy is converted to 
kinetic energy, the resulting kinetic energy E'j^ = Ek + 
Ej can be used to determine when the MI forms. Figure 
|2Ic) shows that E'j^ exhibits a local minimum where the 
Mott core appears in the system. 

An enhanced signal of the features observed in Ej and 
Ej + Ek can be obtained by calculating their numerical 
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derivatives. As shown in Figs. Efb) and[5Ic), a large re- 
duction in the derivatives of both quantities occurs just 
before the MI sets in the center of the trap. The mini- 
mum in Ej + Efc is signaled by a vanishing derivative. 
The derivative of the kinetic energy, Fig. [3a), exhibits 
another signal, a weak maximum and a steady decrease 
after the formation of the MI core. 

Using a Feshbach resonance, one can also obtain a final 
E'j^ exhibiting the features observed in Ej, Fig. ^h). 
To that end, just before turning off the trap, one could 
increase U (and hence Ej) by a large constant factor, 
in a timescale much smaller than the scale set by the 
hopping parameter (h/t). In this case E'^, measured 
after the expansion on the lattice, would be dominated 
by Ei thus allowing the observation of the interaction 
strength at which the MI appears in the trapiS^ 




FIG. 3: (color online). Ground-state on-site interaction en- 
ergy per lattice site (e/) of fermions as a function of the den- 
sity [obtained from a periodic system (no trap) with 102 lat- 
tice sites], for different values oi U/t. At n = 1 fermions are 
in the MI phase, and for any other value of n they are in a 
metallic phase. 

In order for the above proposal to work, one needs a 
fast conversion of the interaction energy into kinetic en- 
ergy. We expect this to be the case since the reduction of 
the density during the expansion produces a preponder- 
ance of empty and singly occupied states over the dou- 
bly occupied ones, i.e., on-site interactions get strongly 
suppressed. For example, in the ground state, the de- 
pendence of the on-site interaction energy on the density 
is shown in Fig. O One can see that below n = 0.25, 
the interaction energy is negligible, independent of the 
on-site repulsion U. 




FIG. 4: (color online), (a) Trap (ETrap) and total (Erotai) 
energies corresponding to the systems of Figs. |2Id)-|2If). (b) 
Ek + Ei for systems with the same trap parameters of Fig. 
|5| and different inverse temperatures l3 = l/ksT. pt = 64 
corresponds to the ground-state results of Figs. |2Id)-|2If). 



Measuring energies to detect the formation of MI do- 
mains has the advantage of being useful independent of 
the statistics of the particles involved. In Figs. |2td)- 
|2If), we also show results corresponding to the bosonic 
system in Fig. The same features discussed above for 
fermions are present in the bosonic case. Wc should stress 
here that the minimum that signals the formation of the 
MI in Fig. I^f) is an exclusive property of trapped sys- 
tems, i.e., it is not present in the periodic case. This 
minimum is related to the fast increase of the trapping 
energy i^xrap [see Fig.0Ja)] produced by the formation of 
Mott domains, which push particles to the outer regions 
of the trap. On the other hand, as expected from the 
absence of a global phase transitionjii the total energy 
of the system increases smoothly throughout the forma- 
tion of the Mott plateau [Fig.^a)]. Hence, we expect the 
above signatures to be present also in higher-dimensional 
systems where most of the experiments are carried out. 

To what extent will finite temperatures reduce this sig- 
nal? We find that as long as the temperature is low 
enough to allow for the formation of a MI phase, the 
dips in the energy are present. In Fig. EJb), we show 
Ek + Ej for systems at different inverse temperatures 
(3 — l/ksT (where is the Boltzmann constant and T 
is the temperature). For fit = 64, the system is essen- 
tially in its ground state [results in Figs. [2Id)-[2If)]. As 
seen in Fig.^b), the ground-state minimum of Ek+Ej is 
still present (although weakening) with increasing tem- 
perature {(3t = 8 and 4) when Mott domains are still 
observed in the density profiles. For high temperatures 
[/3t = 2 in Fig. 01b)] , no Mott plateau appears in the trap 
and no minimum is seen in Ek + Ej . 
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FIG. 5: (color online). Intensity plot of the density fluctua- 
tions Ai = (nf) — (rii)^ vs position and interaction strength 
in fermionic (a) and bosonic (b) systems with 40 atoms on 60 
latice sites and V2a^ — O.Oli. With increasing [/, two Mott 
domains appear at the sides of the central metallic phase be- 
fore a full Mott core develops in the center of the trap. In 
(b), for low values of U /t, Ai has been truncated as in Fig.Q 

In general, Mott-insulating plateaus can also appear 
on the sides of a central metallic region with n > 
With increasing U, these MI regions spread inward and 
merge to occupy the core of the trapiSi2iiLi£ The evo- 
lution of Ai for a system like this is shown in Fig. [3 
for fermions and bosons with the same trap parameters. 
The emergence of two Mott domains around the central 
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FIG. 6: (color online), (a), (d) Kinetic energy (Ek), (b), 
(e) interaction energy (Ei), and (c), (f) the sum Ek + Ei, 
for fermions (a-c) and bosons (d-f) in the systems of Figs. |^ 
(40 atoms and V2a^ = O.Olt). Continuous lines depict the 
numerical derivatives of the energies. The arrows indicate 
when two MI domains (lower value of U/t) and the MI core 
(larger value of U /t) emerge in each system. 



metallic phase is reflected by a suppression of A around 
U/t ^ 3.8 for fermions (a), and U jt ^ 6.3 for bosons (b). 
The central MI core develops at U/t ~ 5.0 for fermions 
(a) and U/t ^ 7.1 for bosons (b). 

In Fig. 13 wc show the kinetic and interaction energies 
and their sum for fermions [(a)-(c)] and bosons [(d)-(f)]. 
Clear local minima can be seen in Ek+Ej, Figs.EIc) and 
ISff), when (i) The two MI domains of fermions (bosons) 
start to develop at the sides of the central metallic (su- 
perfluid) phase, and (ii) when the Mott core forms in 
the center of the trap due to the merging of the two 
Mott shoulders. An increase of Ek + Ej between these 
two points occurs because when the two MI domains sur- 
round the central metallic (superfluid) region, the trans- 
fer of particles to the outer regions costs a finite amount 
of trapping energy. This produces a "freezing" of the 



density profiles with increasing U so that Ek + Ei in- 
creases by finite amounts before particles can be trans- 
fered out of the Mott regions to produce a full MI in the 
trap centeriSi We have also found that the behavior of 
the energies between points (i) and (ii) depends on the 
trap parameters. This means that complicated structures 
may be observed related to the discontinuous transfer— 
of particles out of the central metallic (superfluid) region 
[see for example the fermionic case in Figs.E^b) andlBfc) 
in contrast to the bosonic case in Figs. EJe) andEff)]. 
They are not related to the formation of new Mott do- 
mains. One can be certain that the full Mott plateau in 
the center of the trap has been formed from the smooth 
increase of Ek + Ei with U after a minimum. (In this 
case, the density profiles almost do not change with 

In Fig. El wc have also shown the numerical derivatives 
of the energies. They also provide clear evidence for the 
formation of Mott domains. In the fermionic case, the 
derivative of Ek [Fig. El^a)] decreases steadily after the 
full Mott insulator sets in the middle of the trap. Two 
very weak maxima can be also seen at the points signaled 
by the arrows. [They are clearer in the bosonic case, 
Fig. EJb).] The two values of U/t where the derivative 
of Ek + El, Figs. EJ^c) andElO' increases through zero 
provide a further signature of the two minima in Ek+Ej. 

In summary, we have shown that minima in Ek + Ej 
signal the formation of Mott domains of fermions or 
bosons in a trap. We have also studied the signatures 
that MI phases imprint in the derivatives of the kinetic 
and interaction energies, and of their sum. We have dis- 
cussed the relevance of our findings to the detection, by 
means of TOF measurements, of the formation of Mott 
domains of fermions and bosons confined on optical lat- 
tices. Our proposal represents a novel way to detect 
the metal-Mott-insulator transition using TOF measure- 
ments in fermionic systems, while for bosons the imple- 
mentation of our findings provides an alternative way 
to support the established experimental evidence of the 
supcrfluid-Mott-insulator transition. 
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